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ABSTRACT 



Context. One class of protoplanetary disc models, the X-wind model, predicts strongly subkeplerian orbital gas velocities, a configu- 
ration that can be sustained by magnetic tension. 

Aims. We investigate disc-planet interactions in these subkeplerian discs, focusing on orbital migration for low-mass planets and gap 
formation for high-mass planets. 

Methods. We use linear calculations and nonlinear hydrodynamical simulations to measure the torque and look at gap formation. In 
both cases, the subkeplerian nature of the disc is treated as a fixed external constraint. 

Results. We show that, depending on the degree to which the disc is subkeplerian, the torque on low-mass planets varies between the 
usual Type 1 torque and the one-sided outer Lindblad torque, which is also negative but an order of magnitude stronger In strongly 
subkeplerian discs, corotation effects can be ignored, making migration fast and inward. Gap formation near the planet's orbit is more 
diflrcult in such discs, since there are no resonances close to the planet accommodating angular momentum transport. The location of 
the gap is shifted inwards with respect to the planet, leaving the planet on the outside of a surface density depression. 
Conclusions. Depending on the degree to which a protoplanetary disc is subkeplerian, disc-planet interactions can be very different 
from the usual Keplerian picture, making these discs in general more hazardous for young planets. 

Key words, planets and satellites: formation - planetary systems: protoplanetary discs 



1. Introduction 

Planets form in circumstellar discs, and the gravitational interac- 
tion between planet and disc plays a major role in shaping plan- 
etary syste ms. Tidal waves excited by th e planet lead to orbital 
migration jGoldreich & Tremaine 1979h. which comes in sev- 
eral flavours. Type I migration (' Wardll997h is thought to hold for 
low-mass planets, up to a few times the mass of the Earth (M®). 
High-mass planets, comparable to Jupiter (with mass Mj), can 
\ tidally truncate the disc and open up a deep annular gap around 
i their orbit. The resulting Type II migration is drive n by the vis- 
■ cous evolution of the disc ( Lin & Papaloiz ou 1986). A strongly 
dynamical form of migration. Type III (Masset & Papaloizou 
. 120031) . can be achieved for intermediate-mass planets that par- 
tially open up a gap in mas sive discs, or by pla nets placed 
on a strong density gradient dPepliiiski et al.ll2"008h. For a re- 
cent o verview of disc-planet interactions, see I Papaloizou et alj 
(12001) . 

Most studies of disc-planet interactions have focused on 
discs that have a Keplerian velocity profile, possibly with a slight 
correction for a radial pressure gradient. One class of protoplan- 
etary disc models, however, predicts strongly subkeplerian mo- 
tion d ue to magnetic tensi on jShu et alj|2007l) . In an interesting 
study, lAdams et alj ( l2009l) point out that in these X-wind discs, 
migration due to the strong head wind experienced by the planet, 
which is moving with Keplerian velocity, can overcome Type I 
migration for Earth-sized planets. On the other hand, gap for- 
mation should not be affected that much, with the planet located 
slightly towards the inner edge of the gap. 

In this Letter, we study disc-planet interactions in subkeple- 
rian discs, under the simplifying assumption that the subkeple- 



rian nature of the disc is an external and fixed constraint. This 
allows us to use purely hydrodynamical models as a first ap- 
proximation. Future models should self-consistently include the 
evolution of the magnetic field. We outline the numerical and 
physical set-up in Sect. |2] present the results in Sect.[3]and con- 
clude in Sect.m 



2. Disc model 

The basic equations are conservation of mass, momentum, and 
energy in a 3D cylindrical geometry (r, (p, z). We make two sim- 
plifications at this point: by considering an isothermal equation 
of state, we remove the need to solve the energy equation, and 
by integrating the resulting system vertically, we end up with 
a two-dimensional problem. Although it has been pointed out 
that the isothermal approximati on is not valid in the inner re- 
gions of protoplanetary discs flPaardekooper & Mellemal[2006al 
2008), we see below that in strongly subkeplerian discs these 
effects play essentially no role. The sound speed is given by 
Cs - HQ.(), where H is the scale height of the disc and Qo the 
equilibrium angular velocity (see below). We varied h - H/r^, 
where rp is the orbital radius of the planet, between 0.05 and 
0.2. The higher values of H are appropriate if the disc close to 
the star is 'puffed up' by direct iiTadiation. A kinematic viscos- 
ity V was used, with v(r) such that the initial, constant surface 
density is a stationary solution with no accretion. In the text, 
we quote the corresponding cc-parameter at the location of the 
planet, a - v(r^)/(H^ Q.o). The resulting set of eq uations can 
be found elsewhere (see lPaardekooper & Papaloizou! 2009). The 
value of the surface density at r = rp, Sp, can be chosen arbitrar- 
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Fig. 1. Position of the inner (open triangles) and outer (solid cir- 
cles) Lindblad resonances in discs with varying /. Pressure ef- 
fects have been ignored in calculating the resonance positions. 



ily in non-selfgravitating discs; we normalise the torque accord- 
ingly (see Sect. 13. lb . 

The gravitational potential entering the equations contains 
terms due to the central star and the direct and indirect com- 
ponent of the planet potential. We use a softened point mass 
pote ntial for the planet, with a soften ing parameter b - Q.6H 
(see iPaardekooper & Papaloizoull200 9). This value of b is ap- 
propriate to account for 3D effects. The planet is assumed to be 
on a fixed circular orbit. For the central star potential we use 
0» - -f^GMt/r, where / is a dimensionless parameter govern- 
ing the degree to which the disc is subkeplerian. It is easy to see 
that the resulting equilibrium angular velocity, neglecting pres- 
sure effects, is given by Qo = f^K, where Qk is the Keplerian 
angular velocity. For / - 1, we have a purely Keplerian disc, 
while for / < 1, the disc is subkeplerian. X-wind models sug- 
est / X 0.6, o r even / 0.34 during FU Orionis outbursts 
Shu et al.ll2007h . 

W e combine a linear code jPaardekoope r & Pa paloizoul 
|2009|) wi th fully nonlinear hydrodynami c calculations with 
RODEO (IPaardekooper & Mellernall2006b l) to solve the result- 
ing set of equations. The computational domain extends from 
r/rp - 0.4 to r/r^ - 2.5, and covers the full 2jt in azimuth. This 
domain is covered with a regular grid with 256 cells in the radial 
and 768 cells in the azimuthal direction. Nonreflecting bound- 
ary conditions were used. We vary the planet mass, Mp, and / 
to study different regimes of disc-planet interactions in subkep- 
lerian discs. 



3. Results 

The most important difference between subkeplerian discs and 
Keplerian discs, as far as disc-planet interactions are concerned, 
is the change in positio n of the Lindblad and c orotation reso- 
nances. It was shown in iGoldreich & Tremaind (fT979i) that an 
embedded object can exchange angular momentum with the disc 
only at the location of a resonance. Corotation resonances oc- 
cur where the disc corotates with the planet, and inner (outer) 
Lindblad resonances are located inside (outside) corotation. The 
strong est Lindblad resonances occur at a distance ~ H of the 
planet f^rdl (1 19971) . In Keplerian discs, a planet on a circular or- 
bit is located close to corotation, with Lindblad resonances on ei- 
ther side of its orbit. In strongly subkeplerian discs, the distance 
of the planet to corotation is large enough so that we can safely 
neglect any nonbarotropic effects associated with the corotation 




Fig. 2. Absolute value of the torque on a ^ = 1.26 ■ 10"^ planet 
embedded in discs of constant initial surface density for different 
values of h and /. Solid curves denote the results of linear calcu- 
lations, while the symbols denote results from hydrodynamical 
simulations. 



torque (IPaardekooper & Papaloizoull2008l) . In this sense, subke- 
plerian discs are easier to handle, since we only have to deal with 
Lindblad resonances. 

The approximate position of the Lindblad resonances up to 
10th order are shown in Fig.[T]for different values of /. Pressure 
effects have been ignored in calculating these positions. In a 
purely Keplerian disc if - 1), the planet interacts with both 
inner and outer resonances. For / < 1, the resonance positions 
move inward, with the result that the planet now mainly interacts 
with the outer Lindblad resonances. This has consequences for 
both migration and gap formation. 

3.1. Migration of iow-mass planets 

We first consider the torque F on low-mass planets, not massive 
enough to open up a gap. We compare results from linear calcu- 
lations to hydrodynamical simulations for different values of / 
in Fig.|2] The torque is normalised by Fq = q^l,prpQ.p/h^, where 
q = Mp/M» and Qp is the planet's angular velocity. We plot the 
absolute value of the torque; it is negative for all values of /. 

It is immediately clear from Fig. |2| that the torque very 
strongly depends on /. This can be understood in terms of a 
change in the position of the resonances with respect to the 
planet. For / = 1 we find the class ical Type I torque, modified 
by nonlinear effects at corotation ( IPaardekooper & Papaloizoul 
2009). Around 1 - / h, the planet mainly interacts with the 
strongest outer Lindblad resonances, yielding a very strong neg- 
ative to rque. The to tal torque approaches the one-sided Lindblad 
torque (IWardll997h . which would be obtained when considering 
only outer Lindblad resonances and which is a factor of ~ 1/h 
stronger than the differential Lindblad torque. The maximum 
torque depends on the adopted smoothing, with higher values for 
lower values of b. For 1 - / » /i, the torque decreases, since the 
resonances it interacts with become weaker as the planet moves 
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further away from corotation. For h - 0.05 and / = 0.6, the re- 
sulting torque is comparable in magnitude to the classical Type I 
torque. Linear calculations and nonlinear simulations agree very 
well, which again indicates that nonlinear effects at corotation 
do not play an important role. 

The results depicted in Fig.|2]indicate that the crucial param- 
eter determining the total torque on the planet is p = (1 - f)/h. 
The maxima of the curves occur factors of 2 apart in Fig.|2] con- 
firming this picture. The scaling of the torque with Fq ensures 
that for / = 1, the curves fall on top of each other. When p ~ I, 
a sizeable fraction of the one-sided Lindblad torque can be ex- 
erted on the planet. The one-sided torque scales as h^^, while the 
differential Lindblad torque scales as h^^. This is the reason for 
the different maxima for different values of h. This scaling holds 
for all / < 1 ; only for / = 1 do we get the h^^ scaling of the dif- 
ferential Lindblad torque. Therefore, in subkeplerian discs, for a 
fixed value of p, the torque scales as h^^. For a fixed value of /, 
the situation can be different. For example, for I - f - 0.4, the 
torque for h = 0.2 is twice as strong as would be expected from 
an h^^ scaling, which is due to the shift of the maximum torque 
to higher values of /. 

The torque due to a headwi nd in a subkepleria n disc on a 
planet with radius Rp is given bv I Adams et a (120091) : 

Tx/ro-^Cuil-f)'-J^]\ (1) 
2 q^\rpl 

with Cd a constant of order unity (IWeidenschilUngiri977l) . This 
means that the ratio of the torque due to a headwind and the 
classical Type I torque is proportional to h, and is of order unity 
for an Earth mass planet located at 1 AU, in a disc with / = 0.66 
and h = 0.05. Since most resonances are located far inward with 
respect to the orbit of the planet, the actual Type I torque is close 
to its classical value at / = 1 (see Fig. |2]l. It is however very 
difficult to generalise this result due to the complex behaviour of 
the Type I torque when varying h and /. For example, while for 
/ - 0.6, Fx ~ r at 1 AU for a planet of 1 in a disc with h - 
0.05, the Type I torque will in fact dominate for higher values of 
h. The dimension of the parameter space is quite high (q, rp, h, 
etc.), so that it is difficult to make more general statements. 

It is important to realise that subkeplerian discs are haz- 
ardous environments for low-mass planets. Since the planet 
mainly interacts with outer Lindblad resonances, and corotation 
torques play a minor role, migration will be directed inward. The 
time scale for migration can be even smaller than the classical 
Type I time scale, in unfavourable cases more than an order of 
magnitude smaller 

3.2. Gap formation 

We now turn to the issue of gap formation. The presence of a gap 
determines whether a planet moves in the Type I or the Type II 
regim e, the latter being a much slower mode of migration (WarJ 
Il997h . There are two requirements that have to be fulfilled to 
open a gap: angular momentum transport due to the presence of 
the planet must exceed the viscous transport, otherwise the gap 
could be refilled by viscous evolution of the disc, and, in discs 
with low viscosity, the waves transporting angular momentum 
should be damped locally. For standard parameters, both criteria 
give similar minimum masses for gap o pening, of the ord er of 
10 Me. A unified criter ion is derived in lCrida et~an (l2006h . 

lAdams et al.l (l2009h derived a minimum mass for gap- 
opening in subkeplerian discs, q^^^ - 7t(1 - f)^ah^. They noted 
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Fig. 3. Azimuthally averaged surface density, in units of the ini- 
tial (constant) surface density So, for a ^ = 5 ■ 10"'* planet em- 
bedded in a /z = 0.05, a - 0.001 disc for different values of /. 
The inset shows the averaged surface density at the orbital radius 
of the planet as a function of /. 



that for / = 0.66 this mass is comparable to the classical esti- 
mates for Keplerian discs, although the planet would be located 
more towards the inner edge of the gap. In this section, we study 
gap formation for a ^ = 5 • 10"^ (corresponding to 0.5 Mj around 
a 1 Mq star) planet embedded in a disc with h - 0.05 and a 
viscosity corresponding to a = 0.001. For 0.4 < / < 1, this 
planet fulfils the criterion for gap formation as given above, if 
only slightly for / = 0.4. 

The results are depicted in Fig. [3] where we show the az- 
imuthally averaged surface density for different values of / af- 
ter 500 orbits of the planet. The system has reached a steady 
state by that time. For a Keplerian disc, the surface density close 
to the planet has dropped by more than an order of magnitude, 
confirming that this planet indeed opens up a gap. The main dif- 
ference between Keplerian and subkeplerian discs is the posi- 
tion of the gap, which is shifting inward for / < 1 . This can 
again be understood in terms of resonance positions: the impor- 
tant resonances are located inside the orbit of the planet. The 
associated waves are damped locally through shocks, which re- 
sults in a gap at r < rp. Therefore, contrary to the findings of 
lAdams et al.l (|2009|), the planet is locate d near the outside of the 
gap. The analysis of lAdams et al.l (|2009|) only applies to the case 
1 - f < h, for which the location of the original gap (for / = 1) 
is still densely populated by resonances. In this case, because the 
outer resonances are located closer to the planet (but still outside 
its orbit), the outer part of the gap forms more easily, leaving the 
planet closer to the inner edge of the gap. 

The depth of the gap does not strongly depend on /. As long 
as 1 - / is not larger than the width of the gap for f - I, the 
angular momentum flux will be similar in both cases. Then, the 
ability of a planet to open up a gap is relatively independent of 
/. However, from the perspective of the planet the case / < 1 
can be radically different from the Keplerian case. For what ac- 
cretion and migration are concerned, what matters is the surface 
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Fig. 4, Azimuthally averaged surface density, in units of the ini- 
tial (constant) surface density So, for a q - 10 "* planet embed- 
ded in a /i = 0.05, a - 0.004 disc for different values of /. 



density close to the planet. Since the planet is located more and 
more outside the gap for low values of /, it could accrete and 
migrate as if there were no gap. This is illustrated in the inset of 
Fig. [3] where we show the azimuthally averaged surface density 
at three different locations in the disc for different values of /. In 
Keplerian discs, planets usua lly interact with material within 2 
scale heights from their orbit dBate et alj|2003h . This is the max- 
imum distance that is considered in the inset of Fig. [3] 

For a Keplerian disc, a gap forms with a half width that is 
approximately 2H, and in the inset of Fig. [3] all three curves fall 
below 0.1 for f - 1. For / - 0.8, there is still a density de- 
pression around the orbit of the planet, but for lower values of / 
there is no clear gap near r = rp, it has shifted inwards enough so 
that the planet is basically embedded in the disc again. However, 
since for the lowest values of / there are no more resonances 
located close to the planet, the Type I torque will not be fully 
restored. We find that typically r(/ = 0.6)/r(/ = 1) « 2 for 
gap-opening planets. The influence of the head wind (see Eq.[Tll 
will be strong, however, since it does not rely on resonances. We 
comment that a gap-opening planet is still tidally locked to the 
gap, just as in the Keplerian Type II migration case. 

The possibility of a high-mass planet fully embedded in the 
disc may have some important consequences for gas accretion. 
In a Keplerian disc, a ccretion drops by an o rder of magnitude 
when a gap is formed ( iD'Angelo et alj|2003 l). In a subkeplerian 
disc, there a significant amount of mass remains near the orbit of 
the planet, making accretion potentially very efficient. It is not 
clear, however, if the planet is able to accept material that has 
such a high relative velocity. 

The results for a more massive planet of <7 = 0.001 (1 Mj 
around a 1 Mq star) are very similar to those presented in Fig. [3] 
The gap shifts inward, and for / < 0.7 the planet is located on 
the outer edge of its own density depression. For smaller plan- 
ets, which only open up a shallow density depression for f - I, 
remain fully embedded for / < 1, as shown in Fig.H] While for 
/ = l,aq = 10""* planet decreases the surface density around its 
orbit by a factor 0.7, for / < 0.7 there is no more evidence for 
any density depression. For this lower planet mass, the impor- 
tant resonances become too weak to affect the surface density in 
strongly subkeplerian discs. 

4. Discussion and conclusions 

We have presented hydrodynamical simulations of planets em- 
bedded in subkeplerian discs. They represent the first step 
towards modelling fully 3D magnetised discs. It is known 



for Keplerian discs that magnetic fields can have a strong 
impact on planet migration: regular fields introduce mag- 
netic resonances (Terquem 2003), while magnetic turbulence 
introduces stochastic migration dNelson & Papaloizoul 120041 : 
lAdams & Blochl2009l) . It remains to be seen what impact a mag- 
netic field configuration that gives rise to subkeplerian discs can 
have on the simple hydrodynamic picture presented here. 

We have worked in the isothermal limit, but since corota- 
tion torques only play a minor role in strongly subkeplerian 
discs, results for more realistic discs should be similar The 
two-dimensional approximation, in combination with a grav- 
itational softening parameter of order h, gives similar results 
to fully three-dimensional simulati ons, again at least as far as 
the Lindblad torque is concerned dPaardekooper & Papaloizoul 
2009). We have found that the results do not depend strongly on 
the initial surface density profile. 

We have considered migration of low-mass planets (the Type 
I regime), finding that there is a strong dependence on /. For 1 - 
f ~ h, the planet feels almost the full one-sided Lindblad torque, 
which is a factor l/h stronger than the classical Type I torque. 
Such a disc would be very hazardous to low-mass planets, since 
inward migration is sped up by more than an order of magnitude 
compared to Keplerian discs. For \ - f > h, the torque decreases 
because the resonances the planet interacts with become weaker 
The dependence of the torque on h and / is quite complicated, 
and it is not easy to say for a disc of given / and h whether 
the Type I torque will be stronger of weaker than the head-wind 
torque. 

Gap formation proceeds similar to that in Keplerian discs. 
However, because of the inward shift of the important reso- 
nances, the gap will be located inside the planet's orbit. This then 
leaves the planet on the outside of its own gap. For strongly sub- 
keplerian discs, a gap-opening planet can become fully embed- 
ded again. Since the torque-generating resonances are located 
far away, this does not restore the full Type I torque. Accretion 
time scales could be very short, if the planet is able to accept the 
available matter. 
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